Spin-orbit coupling and magnetic spin states in cylindrical quantum dots 
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We make a detailed analysis of each possible spin-orbit coupling of zincblende narrow-gap cylin- 
drical quantum dots built in a two-dimensional electron gas. These couplings are related to both 
bulk (Dresselhaus) and structure (Rashba) inversion asymmetries. We study the competition be- 
tween electron-electron and spin-orbit interactions on electronic properties of 2-electron quantum 
dots. 
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The creation and manipulation of spin populations in semiconductors has received great attention since the Datta- 
Das proposal of a spin field-effect transistor^ based on Rashba spin-orbit coupling of electrons in a bidimensional 
electron gasj^ and the possibility for quantum computation devices using quantum dots (QDs). 3 Thus, it is important 
that every spin-orbit (SO) effect be clearly understood for a full control of spin-flip mechanisms in nanostructures. 

There are two main SO contributions in zincblende materials. In addition to the structure inversion asymmetry 
(SIA) caused by the 2D confinement (the Rashba SO), there is a bulk inversion asymmetry (BIA) term in those 
structures (the Dresselhaus SO)^ An yet additional lateral confinement defining a dot introduces another SIA term 
with important consequences, as we will see in detail. Although the relative importance of these two effects depends 
on the material and on sample design (via interfacial fields), only recently have authors begun to consider the behavior 
of spins under the influence of all effects. 

The goal of this work is to show how important different types of SO couplings are on the spectra of parabolic 
QDs built with narrow-gap zincblende materials. We consider both Rashba and a diagonal SIA, as well as the all 
Dresselhaus BIA terms in the Hamiltonian, in order to study features of the spectrum as function of magnetic field, 
dot size, and electron-electron interaction. 

Consider a heterojunction or quantum well confinement potential V(z) such that only the lowest z-subband is 
occupied. The Hamiltonian for a cylindrical QD, in the absence of SO interactions, is given by Ho = (fi, 2 /2m)k 2 + 
V(p) + 5^ B B • er/2, where k = — iV + eA/(hc), A = Bp(— sin 9, cos 6, 0)/2 describes a magnetic field B = Bz, m 
is the effective mass in the conduction band, g is the bulk <?-factor, p B is Bohr's magneton, V(p) — mur^p 2 /2 is the 
lateral dot confinement, and cr is the Pauli spin vector. The analytical solution of Hq yields the Fock-Darwin (FD) 
spectrum, E n i az — (2n + \l\ + l)fi£l + lhwc/2 + gp, B Boz /2, with effective (cyclotron) frequency = \/ujq + u 2 c /4 
(u>c = eB / (mc)). The FD states are given in terms of Laguerre polynomials. 5 ! The lateral, magnetic and effective 
lengths are l = y/hj (muj ), l B — \/K/ {mu>c) and A = ^/h/(mil), respectively. 

The SIA terms— for the full confining potential, V(v) = V(p) + V(z), and coupling parameter a are decomposed as 
Hsia = Hr + Hg IA + Hk- These three forms are: i) Hk = ia(htuo/lo)Xx[(T + L- — <j-L + ] (k z ) gives zero contribution 
when (k z ) = (pure state parity); ii) Hg IA = aa z(hu>o /Iq)[L z + X 2 x 2 /2l B ] is the diagonal contribution due to 
the lateral confinement and x — p/X been a dimensionless radial coordinate, Lz — —id/ 88 is the z-orbital angular 
momentum; iii) Hr — —a(dV/ \dz)[a+L-A- + a-L+A+] is the Rashba term for the perpendicular confinement 
dV/dz, L± — e ±l9 , a± = (ax ± iay)^ and A± = \^d/dx + Lz/x + x\ 2 /(2l%)]. In principle these terms can be 
tuned since Hg IA depends on the confining frequency wo while Hr depends on the interfacial field dV/dz. 

The BIA Hamiltonian 4 for zincblende materials, after averaging along the z-direction, is given by Hbia — 
(k 2 ) [o-yky - a x k x ] + 7cr 2 (k z ) (k 2 - k 2 ), where 7 is the Dresselhaus parameter. The first 
(second) term is cubic (linear) in the in-plane momentum. The last term will be zero for systems where (k z ) = 0, 
while (k 2 ) — (tt/zq) 2 , zq being the z-direction (perpendicular) confinement length. In cylindrical coordinates H B ja 
can be written as H B ja = H^ + H£, where H^ = — i^/X) [a + L + A + — <j-L-AJ\ (fc 2 ) is the linear term and, after 
long algebra manipulation, 6 the cubic term becomes H^ — i(j/X 3 ) \a-L\H\ + a + L i _H2 + (7-L-H3 + a+L+H^, 
where Hi = Ai + jrBi + jrCi + j$-Di, with i = 1,2, 3, 4. The long expressions for the sixteen functions Ai, Bi, Ci, 

— — ''B B 

Di are given in Ref. |3- 

Finally the electron-electron interaction H ee ~ e 2 /[e|ri — 1*2 1] , with e being the dielectric constant of the material, 
is expanded into Bessel functions Jfc(£) as H ee — (MlX/as) Sfcl-00 e lfc ^ 1 ~ 6 ' 2 ^ J °° d£ Jk(£,xi) Jk(£,X2)e~^ z °t x , where 
flB = eh 2 1 (me 2 ) is the effective Bohr radius. The FD basis states must be properly antisymmetrized to describe 
unperturbed spin eigenstates£ 

Summarizing, our total single-particle Hamiltonian is given by H = Hq + Hg IA + Hr + Hp + H£. For the 
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two-particle case, we study the states and spectrum of H + H ee . Parameters for InSb are in Ref. [7j. 

We present results by analyzing the role of each SO term in the Hamiltonian. We take into account all states in the 
FD basis having n < 4 and \l\ < 9 in our numerical diagonalization, which is equivalent to the first ten energy shells 
at zero field and embodies a total of 110 basis states. The sequence of FD states of Hq starts at zero B- field with 
{n,l,(Tz} = {0,0,±1}, followed by the degenerate set {0, —1, ±1} and {0, 1, ±1}. The next energy shell is composed 
by {0, —2, ±1}, {1, 0, ±1}, and {0, 2, ±1}£ Spin and orbital degeneracies are broken by B and the states with negative 
I and positive oz acquire lower energies because of the negative g-factor of InSb. The lowest FD crossing occurs 
between states {0, 0, —1} and {0, —1, +1}, at a critical field B c — mhw /[ii B ^m\g\{fh\g\ + 2)]. The moderate value 
of B c in InSb (~ 2.6 T for Hujq = 15 meV) is a direct consequence of its large \g\ factori For GaAs = 0.44, 
m = m/rriQ = 0.067), this level crossing would appear only at B c (GaAs) ~ 9.5 T for a much smaller confinement, 
Hujq = 2 meV, corresponding to a regime where Landau levels are well defined. 

Any figure showing spectrum has the structure: Panel A shows QD spectrum for the full FD basis (110 states); 
Panel B shows a zoom on the three lowest shells, plus inset with another zoom on the 4 levels of the second shell; 
Panel C and D show, respectively, the devolution of spin oz and orbital I angular momenta, for the full FD basis, 
while their insets take into account only the lowest 7 QD levels. 

Figure 1 shows the simultaneous addition of both SIA terms, Hg IA and Hr, to Hq. The diagonal term Hg IA causes 
small splittings on the zero-field spectrum and readjust the sequence of states according to total angular momentum 
j = l + az/2. For example, the highest (lowest) zero-field energy level in the second shell has j — 3/2 (j — 1/2). Since 
Hg IA does not induce shift on the accidental degeneracy points of the FD spectrum at finite fields, the first level 
crossing occurs at a critical value Be — B c ~ 2.6 T. Also, this term does not induce any level mixture on the FD 
states. The Rashba term Hr introduces a strong state mixture for any magnitude of a. This is evident for any pair 
of FD levels satisfying Al = — Acr^/2 = ±1 that show a crossing at given accidental degeneracy of the FD spectrum. 
The induced mixture converts this crossing at B c to an anticrossing (AC), at a shifted critical field Be — 2.5 T < B c , 
with an energy minigap. Higher energy levels, also satisfying this selection rule, present ACs around the same value 
of Be, and gives origin to the observed collapse in both az and I quantum numbers at B ~ 2.5 T, shown in Panels C 
and D. The range of critical fields (between 2.1 and 2.6 T), and the size of the minigaps opened at those AC regions, 
are proportional to the magnitude of a. Hr also induces small splittings in the zero-field spectrum and slightly shifts 
the accidental degeneracy points at finite fields. After adding both SIA terms simultaneously (full spectrum in Panel 
A), one can see in the inset of Panel B that the ordering of states is the one determined by Hg IA . However, the 
energies of states j = 3/2 at 30 meV and j = 1/2 at a slightly smaller value, as well as the value of field (B > 0.1 
T) where the normal state ordering (the one existing in the absence of SO terms) for g < material is restored, are 
determined by Hr. For increasing energy, the ordering for this second shell is {0, — 1, +1}, {0,-1,-1}, {0,1,+1}, 
{0, 1, —1}. The width of critical fields becomes wider, between 2.2 and 3.6 T, as seen in Panel C. Furthermore, Panel 
D shows that orbitals having I < (I > 0) present ACs at fields smaller (larger) than the field Be — 2.55 T where 
occurred the first AC (see insets of Panels C and D). For a future comparison notice, in Panel B, that the first AC 
near 50 meV involving n = 1 states {0, 1, —1} and {1, 0, +1} occurs at that same Be value. In general, ACs between 
states with any n value occur inside the same unique range of critical fields, as shown in Panels C and D. Finally, 
observe that both SIA terms can be reduced by decreasing uq {Hg IA ) or dV/dz (Hr), and that all even and negative 
I states (I = —2, —4, —6, —8) show anticrossings. 

Figure 2 shows the simultaneous addition of both BIA terms, H^ and Hp, to Hq. The cubic term H£, under 
the present QD parameters, has small influence over the Hq spectrum. Small state mixtures is induced at So — 1 
and ~ 5 T, both involving ACs satisfying A^ = =f3 and Aaz/2 = ±1- The first one, at 1 T (5 T), occurs between 
states {0, 1, —1} and {0, —2, +1} ({0, 0, —1} and {0, —3, +1})- This term also induces zero-field splittings on the FD 
spectrum and also a shift due to matrix elements for Al = ±1 = Aaz/2- However, such splittings and opened minigaps 
at the ACs regions are very small. Therefore, the simultaneous addition of both BIA terms (full spectrum in Panel 
A), where the linear term Hp is the most important, drastically changes the general features of the FD spectrum, 
inducing strong zero-field splittings and shifting its accidental degeneracies to higher fields. Yet, the respective matrix 
elements (Al = ±1 = Act^/2) does not introduce ACs on the lowest energy levels. As seen in Panels C and D, the 
mixing induced by the linear term is so strong that the QD states arc not anymore pure states even at zero field. 
Notice in Panel C that at B = values of \az\ < 0.5 are found for high energy states, while in its inset one finds 
(7 z — 0.7 for the ground state. As an example of level crossings displaced to higher fields, observe in Panel B that 
the first one has moved to Be — 3.3 T, there is only one crossing present in the second shell at about 0.45 T (see 
inset), and the second one occurs a higher field around 3.5 T. Thus, contrary to the observed for SIA case, the normal 
ordering of state is no longer restored. We will come back to this fact later. As a final note, in the same inset and 
at zero field, the highest (lowest) energy state has j = 3/2 with eigenvalue equal to energy of 30 meV of the pure Hq 
(j = 1/2 at smaller energy near 27 meV). The influence of H^ on the spectrum changes with zq. 

Figure 3 shows the one-particle QD spectrum for full H or when all SO terms are simultaneously taken into 
account. From the previous discussions, one may identify which SO mechanism is dominant in each of the main 
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FIG. 1: Spectrum when Hl? IA and Hr are added to Ho (A and B). Critical field range for ACs is seen on C and D. The 
lowest one (Be — 2.55 T, insets) occurs near B c . ACs involving I < (I > 0) orbitals are shifted to lower (higher) fields (D). 



FIG. 2: Spectrum when Hp and Hp are added to Ho- The linear contribution dominates the cubic one. Strong mixing at low 
fields are due to Hp, while the AC around 6 T is due to H£. 

signatures present on the spectrum. An enormous state mixture, even at small magnetic fields (Panels C and D and 
their insets), as also splittings, position and ordering of states (Panels A, B and its inset) are dominated by Hp, 
although with contributions from both SIA terms. The small influence of H£ remains around 6 T. The lowest ACs 
are induced by the Rashba term Hr, although shifted to higher critical fields by the linear BIA term, H^. Observe 
that the first AC has moved from 2.55 T (Fig. 1) to 3.3 T (Panel B and insets of Panels C and D ), and the ranges 
of critical fields becomes wider. However, the new feature of the full H spectrum is the clear presence of more than 
one unique range of critical fields where ACs occur (compare with Fig. 1). Details on Panels A and C: i) The first 
family of ACs, near 3.3 T (related to states between 20 and 70 meV) involves only n — levels. The first ACs being 
between {0, 0, -1} and {0, -1, +1}, {0, -1, -1} and {0, -2, +1}, {0, -2, -1} and {0, -3, +1}, ... ; ii) A second family 
of ACs around 5 T ( related to states between 70 and 120 meV) involves only n = 1 levels. The first ACs being 
between {0,1,-1} and {1,0, +1}, {1,0,-1} and {1,-1, +1}, {1,-1,-1} and {l,-2,+l}, ... ; hi) A third family 
of ACs around 8 T (related to states between 130 and 180 meV) involves only n — 2 levels. The first ACs being 
between {0,2,-1} and {1,1, +1}, {1,1,-1} and {2,0,+l}, {2,0,-1} and {2,-1, +1}, ... . Although lowest ACs in 
the QD spectrum are caused by the selection rules of Hr, the presence of and H£, in the full H, displaces and 
regroup all states with same n value that contribute to the minigap near a fixed critical field value. 

In Fig. 4 we simulate the cancellation of zero-field splittings even in the presence of all SO terms, what is reasonably 
obtained by taking an interfacial field dV/dz four times stronger than that one considered before (see Ref.0, other 
parameters remained unchanged). This is equivalent to increase the influence of the Rashba term Hr. Notice, 



4 



FIG. 3: Spectrum of full H, where Hr induces minigap regions that are shifted to higher fields by Hp. H% has small influence 
but induces state mixtures. The zero-field splittings produced by H® IA are dominated by those from Hp. 

in Panels A and B, that not only the zero-field splittings are nearly vanished, but also the Zeeman splittings are 
practically suppressed at low fields (B < 1.5 T). At zero magnetic field, an energy shell structure identical to the 
pure Hq and with the same level separation of 15 meV is formed at displaced energies. In the inset of Panel B one 
sees that the energy of j — 3/2 level is pushed near j — 1/2 level, going from 30 (in Fig. 3 B) to 26.5 meV (in Fig. 4 B). 
While the zero-field splittings nearly vanish the energy minigaps increases, as seen in Panel B. The rearrangement 
of electronic levels is so remarkable that ACs related to the cubic BIA term at 1.2 T become visible (Panel B and 
insets of Panels C and D). The minigaps at 33 (44) meV involves states {0, 1, —1} and {0, —2, +1} ({1,0, —1} and 
{0, —3, +1}). Even though the electronic levels are less disperse here than in Fig. 3, the SO-induced state mixture is 
much more intense, as can be seen in Panels C and D. Between and 4 T, most of the QD levels have \<rz\ < 0.5 and 
only the ground state has uz — 0.7. As mentioned before, the insets of Panels C and D show that a strong Rashba 
term enlarges the spin- flip region near Be- 



FIG. 4: Full H spectrum with four times stronger dV/dz. Notice the cancellation of zero-field and Zeeman splittings at low 
fields (A, B and inset). New ACs due to H% selection rules occur near 1.2 T (B and insets in C and D). The lowest AC is 
shifted back to Be = 2.7 T. Notice the enormous state mixture in C and, at zero field, most states are displaying \oz\ < 0.5. 

One can further appreciate the intricate balance between SO terms by analyzing how some quantities are affected 
by changes on the lateral and vertical sizes, Iq and zo, or on Rashba field, dV/dz.as shown in Fig. 5. Curves with 
squares, circles and triangles refer to a QD, respectively, with parameters of Ref. 0, with zq doubled (smaller linear 
BIA contribution) and with four times stronger dV/dz, while the dotted curve, in middle Panel, indicates the B^ 
field where the first FD level crossing occurs. The zero-field splitting (left Panel) for states j = 3/2 and j = 1/2 of 
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the second shell is dominated by the linear BIA contribution for any value of Iq- An increase on zq strongly reduces 
the splittings because the Dresselhaus contribution becomes weaker. The reduction is even more drastic by increasing 
dV/dz, which makes Hr larger and, thus, may cancel or suppress zero-field splittings produced by Hg IA . 



FIG. 5: Zero-field energy splittings for the states in the second energy shell (left Panel), critical magnetic fields where the first 
level AC (middle Panel) occurs, and energy minigaps opened at that AC (right Panel) as function of the QD lateral radius 
lo- Meaning of square, circle and triangle curves are explained in text. Arrows at lo = 190 A show the QD radius where the 
spectra from Figs. 1 to 4 were calculated. 

The critical fields Be, where ACs determined by Hr occur (middle Panel, for the lowest minigap between levels 
{0, 0, —1} and {0, — 1, +1}), decrease with increasing QD size, once B c ~ u Q ~ Its value is close to B c when 

BIA terms are not present and the inclusion of Hp shifts Be to higher value. Increasing zq or dV/dz decreases Be, 
once they will decrease the effects due to Hbia- The value Be — 2.1 T, for l a = 270 A (huj — 7.5 meV) is displaced 
to 1.8 (1.5 T) if dV/dz (zq) is four times larger (doubled). This last value can be compared to reported Be — 1.7 
T in Ref. 0, where BIA terms are absent. The small difference AB = 0.2 T can be attributed to the inclusion of 
non-parabolicity effects. Anticrossings at such low fields may be interesting for applications due to easier access. 

Finally, the minigap opened at Be (right Panel) has their main origin in the Hr term, while the inclusion of 
Hbia causes a substantial reduction. If the value of zq is doubled the splitting is enhanced slightly. A yet larger zq 
produces no significant changes. However, the splitting can be drastically enhanced by increasing the Rashba field as, 
for example, changing from 1 to 4.2 meV at Zq = 190 A, when interfacial field is increased four times. Measurement 
of those three quantities should yield important information on the relative strength of SO parameters a and 7. 

After having studied the one-particle QD problem we show, in Fig. 6, the two-electron QD spectrum under 
magnetic field (parameters in 0). On the construction of Slater determinant for two-particle states we use the 20 
lowest one-particle orbitals (|Z| < 3 and n < 1), which amounts to 190 possible two-particle states that can be labeled, 
in the absence of SO interactions, by the projections of orbital (Ml) and spin (Ms) total angular momenta. If no 
SO term is present in H (Panel A), we verified that the singlet ground state is located at 35 meV, while the first 
excited shell at zero-field is splitted by the exchange interaction, being composed by a triplet (at 47.5 meV) and a 
singlet (at 50 meV) states. At very small magnetic field (~ 0.1 T), the normal sequence of QD states is restored. For 
increasing energy and using the notation {Ml, Ms}, the ordering is: {0,0} for the ground state, {—1,1}, {—1,0} , 
{—1,-1}, {1,1}, {1,0} , {1,-1} for the first excited triplet (T), and {—1,0} , {1,0} for the first excited singlet 

(S). The crossing between ground singlet and first excited triplet states occurs at B c ^ =2.1 T. 

Panel B shows the QD spectrum for full two-particle Hamiltonian, H ee + H. We may identify some similar features 
to the single-particle case. For example, the linear Dresselhaus term almost destroys the energy shell structure at 
zero field by shifting level crossings and inducing new zero- field splittings, while the Rashba term introduces energy 
minigaps in the spectrum. Panel B shows details on the competition between Coulomb and spin-orbit interactions in 
narrow-gap cylindrical dots. Observe that the SO interaction, at zero field, acts against the direct Coulomb energy 
and, in a sense, favoring the exchange term. For example: i) The ground state is shifted back from 35 to 31 meV, 
which is close to the non-interacting value 30 meV; ii) The first excited shell states have energies from 43 to 47 meV, 
values even smaller than the non-interacting energy 45 meV. Other important feature in the first excited shell is the 
observation that the original triplet is broken on its three possible terms according to the projection of the total 
angular momentum, Mj — Ml + Ms. For increasing energy order, these terms are composed, at zero field, by the 
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FIG. 6: Two-particle QD spectrum without (A) and with (B) all SO terms. It is seen that the SO energy acts against the 
direct and in favor of the exchange Coulomb energies. The first excited triplet, at zero field, is splitted according to the possible 
Mj values as explained in text. C (Ms) and D (Ml) show the lowest ACs as induced by Hr and shifted by Hp. 



states {-1,1} and {1,-1} (Mj = 0), {-1,0} T and {1,0} T (\Mj\ = 1), {-1,-1} and {1,1} (\Mj\ = 2), in increasing 

s s 

energy, while the ground ({0, 0}, Mj — 0) and first excited ({—1, 0} and {1, 0} , \Mj\ = 1) singlets remain the same. 
Panels C and D show the SO-induced mixing of those lowest states. The first AC at Bq^ = 2.7 T involves states 

{0,0} and { — 1, 1} (Hr selection rule yields AM L = ±1 = — AM$), so that the difference between B^^ and B^ 2e ^ 
(~ 0.6 T) is basically the same between Be and B c (<~ 0.7 T) for the one-electron problem. This means that the 
shifting of Be due to H^ is not altered by QD occupation, although Be itself is decreased by increased occupation. 
Yet, the H^ selection rule becomes AM L = =f3 and AM$ = ±1, and its first minigap between states {—2, 1} and 
{1, 0} S / T would be visible if a higher dV/dz had been considered on the solution for the two-electron problem. 

A very special difference between the one- and two-particle problems is that a strong intrinsic (no phonon-assisted) 
singlet-triplet transition (qubit) at low magnetic fields involving the ground state becomes possible in the two-electron 
case and, in principle, could be explored in implementations of quantum computing devices. As mentioned, the 
critical field is decreased by QD occupation (from Be = 3.3 to B^ 1 = 2.7 T), and this reduction may be increased 
by decreasing the QD confinement energy. At these critical fields where the intrinsic state mixture is enhanced, the 
SO-induced spin relaxation rate (r) can be estimated from the minigap energy (A), as T — h/A. For the lowest AC, 
A values are taken from the right Panel of Fig. 5, from where one sees that A is completely changeable by the QD 
parameters and, consequently, the intrinsic rate T can be changed according to those parameters. 

We showed that inclusion of all SO terms is essential in order to obtain a complete picture of the electronic 
structure of narrow-gap QDs, and discussed the role played by each BIA and SIA terms on QD spectra and on spin 
polarization of states. The combination of strong SO coupling in Hr and large (/-factor introduces strong intrinsic 
mixtures and low excitations on the single-particle spectrum; the position of critical fields where minigaps occur is 
affected by H^. We observed that the two-particle spectrum exhibits strong singlet-triplet coupling involving QD 
ground state at moderate fields, which may have significant consequences like possible use in qubits designs. 
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